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We study two-dimensional (2D) binary-fluid turbulence by carrying out an extensive direct nu¬ 
merical simulation (DNS) of the forced, statistically steady turbulence in the coupled Cahn-Hilliard 
and Navier-Stokes equations. In the absence of any coupling, we choose parameters that lead (a) to 
spinodal decomposition and domain growth, which is characterized by the spatiotemporal evolution 
of the Cahn-Hilliard order parameter (f, and (b) the formation of an inverse-energy-cascade regime 
in the energy spectrum E{k), in which energy cascades towards wave numbers k that are smaller 
than the energy-injection scale kinj in the turbulent fluid. We show that the Cahn-Hilliard-Navier- 
Stokes coupling leads to an arrest of phase separation at a length scale Lc, which we evaluate from 
S{k), the spectrum of the fluctuations of (j). We demonstrate that (a) ~ Lh, the Hinze scale that 

follows from balancing inertial and interfacial-tension forces, and (b) Lc is independent, within error 
bars, of the diffusivity D. We elucidate how this coupling modifies E{k) by blocking the inverse 
energy cascade at a wavenumber fcc, which we show is ~ 27r/Lc. We compare our work with earlier 
studies of this problem. 

PACS numbers: 47.27.-i,64.75.-g,81.30.-t 


Two-dimensional (2D) fluid turbulence, which is of 
central importance in a variety of oceanographic and at¬ 
mospheric flows, is fundamentally different from three- 
dimensional (3D) fluid turbulence as noted in the pio- 
neeri^ studies of Fjprtoft, Kraichnan, Leith, and Batch¬ 
elor [iHg- In particular, the fluid-energy spectrum in 
2D turbulence shows (a) a forward cascade of enstrophy 
(or the mean-square vorticity), from the energy-injection 
wave number kinj to larger wave numbers, and (b) an 
inverse eascade of energy to wave numbers smaller than 
kinj- We elucidate the arrest of phase separation in a 2D, 
symmetric, binary-fluid mixture by turbulence. 


In the absence of turbulence, binary-fluid mixtures 
have played a pivotal role in the development of the 
understanding of (a) equilibrium critical phenomena at 
the consolute point, above which the two fluids mix 
§|, (b) of nucleation [^, and (c) spinodal decomposi¬ 
tion, the process by which a binary-fluid mixture, be¬ 
low the consolute point and below the spinodal curve, 
separates into the two, constituent liquid phases until, 
in equilibrium, a s ing le interface separates the two coex¬ 
isting phases [^, [IJ. In the late stages of growth, as 
the binary-fluid mixture evolves via spinodal decomposi¬ 
tion towards the completely phase-separated, equilibrium 
state, the domains of these two phases coarsen to yield 
ever larger domains whose linear size diverges as a power 
of the time t; this divergence leads to universal sc aling 
forms for the time-dependent correlation functions 


21| of the order parameter (f, which distinguishes the two 


phases of the binary-fluid mixture. 


Coarsening arrest by 2D turbulence has been studied in 
Ref. where it has been shown that, for length scales 
smaller than the energy-injection scale iinj = 2'K/kinj, 


the typical linear size of domains is controlled by the 
average shear across the domain. However, the nature 
of coarsening arrest, for scales larger than linj, i.e., in 
the inverse-cascade regime, still remains elusive. In par¬ 
ticular, it is not clear what happens to the inverse en¬ 
ergy transfer, in a 2D binary-liquid, turbulent mixture, 
in which the mean size of domains provides an additional, 
important length scale. We resolve these two issues in our 
study. By combining theoretical arguments with exten¬ 
sive direct numerical simulations (DNSs) we show that 
the Hinze length scale Lh (see Refs. |23|, [2^) provides a 
natural estimate for the arrest scale; and the inverse flux 
of energy also stops at a wave-number scale ~ 27r/L^f. 
In particular, we study two-dimensional (2D) binary-fluid 
turbulence by carrying out a direct numerical simulation 
(DNS) of the forced, statistically steady turbulence in the 
coupled Cahn-Hilliard and Navier-Stokes equations. In 
the absence of any coupling, our choice of forcing leads 
(a) to spinodal decomposition and domain growth, which 
we examine by the spatiotemporal evolution of 4>, and (b) 
to the formation of an inverse-energy-cascade regime in 
the energy spectrum E(k), in which energy cascades to¬ 
wards wave numbers k that are smaller than the energy- 
injection scale kinj in the turbulent fluid. We show that 
the Cahn-Hilliard-Navier-Stokes coupling leads to an ar¬ 
rest of phase separation at a length scale Lc, which we 
evaluate from S{k), the spectrum of the fluctuations of 
(j). We demonstrate (a) Lc Lh and (b) that Lc is in¬ 
dependent, within error bars, of the diffusivity D. We 
elucidate how this coupling modifies E{k) by blocking 
the inverse energy cascade at a wavenumber kc, which 
we show is ~ 2tt f L^- 

We model a symmetric binary-fluid mixture by using 
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the incompressible Navier-Stokes equations coupled to 
the Cahn-Hilliard or Model-H equations 25, 2^. We are 
interested in 2D incompressible fluids, so we use the fol¬ 
lowing stream-function-vorticity formulation 27-^ for 
the momentum equation: 


idt+u-W)uj = ,^V^uj-W x{cl>Wti) + U, ( 1 ) 

{dt + u ■ V)(j) = and V • it = 0. (2) 


Here u = {ux^Uy) is the fluid velocity, w = (V x 
tt)e 2 , <j){x^t) G [—1,1] is the Cahn-Hilliard order pa¬ 
rameter at the point x and time t, p is the pressure, 
fi{x, t) = 5F[(j)]/5(j){x, t) is the chemical potential, F[(j)] = 
^~ + |V^p/2](ia; is the free energy, A is 

the mixing energy density, f controls the width of the 
interface between the two phases of the binary-fluid mix¬ 
ture, V is the kinematic viscosity, the surface tension a = 
2-\/2/3(A/^), the mobility of the binary-fluid mixture is 
M, and fi^ is the external driving force. For simplicity, we 
study mixtures in which M is independent of (j) and both 
components have the same density and viscosity [l^ . We 
use periodic boundary conditions in our square simula¬ 
tion domain, with each side of length L = 2tt. To obtain 
a substantial inverse-cascade regime, we stir the fluid at 
an intermediate length scale by forcing in Fourier space 
in a spherical shell with wave-number k^nj = ‘2'^/iinj- 
Our choice of forcing fcjiK t) = w(k, t)/Y,k=ki„ 
where the caret indicates a spatial Fourier transform, en¬ 
sures that there is a constant enstrophy-injection rate. 
In all our studies we use kinj = 40 so that there is a 
clear separation between iinj and We conduct DNSs 
of Eqs. dH) and (O by using a pseudospectral method 
because of the cubic nonlinearity in the chemical poten¬ 
tial /i, we use A"/2-dealiasing. For time integration we use 
the exponential Adams-Bashforth method ETD2 (sij. 
Important nondimensional numbers for the turbulent 
flows here are the Grashof number Gr = {L'^ the 

injection-scale Reynolds number Re = Uinjiinj/i^, with 
Uinj = , where tinj is the energy-injection 

rate, the Weber number We = pu^^^iinj/(J^ the Cahn 
number Ch = ^/L, the Peclet number Pe = UrmsL/D, 
where Urms is the root-mean square velocity, and the 
Schmidt number Sc = vjD^ where D = MA/^'^ is the 
diffusivity of our binary-fluid mixture. We give the pa¬ 
rameters for our simulations in the Supplemental Mate¬ 
rial ( 4 ^ . 

Given u{x,t) and from our DNS, 

we calculate the energy and order-parameter 
(or phase-field) spectra, which are, respec¬ 
tively, E{k) = Efe-i<fe'<fc+i(l“(k',i)P)t and 

'S'(fc) = Efe-i<fe'<fc+i(l<^(k':OP)t> where ()t denotes the 
average over time in the statistically steady state of our 
system. The total kinetic energy is E{t) = ^{\\i{x,t)\'^)x 
and the total enstrophy e{t) = 5 (|w(x,t)p)x, where Ox 
denotes the average over space, (/cjw) is the enstrophy- 
injection rate, which is related to the energy-injection 


rate via = {fojU!)/kf^j, E = QA'^j.E{k) is the fluid 
kinetic energy, is the fluid-energy 

dissipation rate, and /c^(|/i(k, t)p|)t is the 

energy-dissipation rate because of the phase field (p. 

Forced, 2D, statistically steady, Navier-Stokes-fluid 
turbulence displays a forward cascade of enstrophy, from 
£inj to smaller length scales, and an inverse cascade 
of energy to length scales smaller than imj- In the 
inverse-cascade regime, on which we concentrate here, 
E{k) ~ (see, e.g.. Refs. ii) and the energy 

flux n(fc) ^ e = assumes a constant value. For 

the Cahn-Hilliard model, if it is not coupled to the 
Navier-Stokes equation, S{k,t) ^ 5(fcL(t)), for large 
times, where the time-dependent length scale L(t) ^ 
in the early Lifshitz-Slyozov [l^ . [isl . [l9l . regime; if 


the Cahn-Hilliard model is coupled to the Navier-Stokes 
equation, then, in the absence of forcing, L(t) ~ t, in 
the viscous-hyd rody namic regime, first discussed by Sig- 
gia 14, 15, [ll, [ 2 J, and L(t) ~ in the very-late- 
stages in the Furukawa [l^ and Kendon 171 regimes. For 
a discussion of these regimes and a detailed exploration 
of a universal scaling form for L(t) in 3D we refer the 
reader to Ref. [T^. We now elucidate how these scaling 
forms for E{k) and S{k,t) are modified when we study 
forced 2D turbulence, in the inverse-cascade regime in 
the coupled Cahn-Hilliard-Navier-Stokes equations. 

In Fig. d] we show pseudo-gray-scale plots of p, at late 
times when coarsening arrest has occurred, for four differ¬ 
ent values of We at Re = 124; we find that the larger the 
value of We the smaller is the linear size that can be asso¬ 
ciated with domains; this size is determined by the com¬ 
petition between turbulence-shear and interfacial-tension 
forces. This qualitative effect has also been observed in 
earlier studies of 2D and 3D turbulence of symmetric 
binary-fluid mixtures 2^, 32-39j|. 



FIG. 1. Pseudo-gray-scale plots of the order parameter field 
(j), at late times when coarsening arrest has occured, in 2D 
symmetric-binary-fluid turbulence with Re — 124. Note that 
the domain size decreases as we increase the Weber number 
We from the leftmost to the rightmost panel: We = 1.2-10“^ 
(R3); We = 5.9 • 10"^ (R4); We = 1.2 • 10”^ (R5); and We = 
5.9- 10”^ (R8). 


We calculate the coarsening-arrest length scale 

L, = 27:[J2S{k)]/[J2kS{k)]. (3) 

k k 

We now show that Lc is determined by the Hinze scale 
Ljj, which we obtain, as in Hinze’s pioneering study of 
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droplet break-up [^, by balancing the surface tension 
with the inertia as follows: 


Lh 


^inj 


( 4 ) 


We obtain for 2D, binary-fluid turbulence the intuitively 
appealing result Lc ^ Lh (for a similar, recent Lattice- 
Boltzmann study in 3D see Ref. [1^). In particular, if we 
determine Lc from Eq. with S{k) from our DNS, we 
obtain the red points in Fig. O which is a log-log plot of 
aLc versus the black line is the Hinze result o 

for Lh^ with a constant of proportionality that we find is 
~ 1.6 from a fit to our data. We see from Fig. [2]that the 
Hinze length scale Lh gives an excellent approximation 
to the arrest scale Lc over several orders of magnitude 
on both vertical and horizontal axes. Note that the Hinze 
estimate also predicts that, for fixed values of and 
CT, the coarsening-arrest scale is independent of D] the 
plot of Lc versus D, in the inset of Fig. [2l shows that our 
data for Lc are consistent (within error bars) with this 
prediction. 

In Fig. [2] (b) we show clearly how the arrest of coars¬ 
ening manifests itself as a suppression of Sik), at small 
k (large length scales). This suppression increases as We 
increases (i.e., a decreases); and S{k) develops a broad 
and gentle maximum whose peak moves out to large val¬ 
ues oi k as We grows. These changes in S{k) are associ¬ 
ated with W e-dependent modihcations in the probability 
distribution function (PDF) P{(j)) of the order parame¬ 
ter which is symmetrical about ^ = 0 and has two 
peaks at ^ = (j)±, where </>+ = —if- > 0 ; we display 
P{cj))/Pjn{cj)) in Fig. [2] (c) in the vicinity of the peak at 
(/)+; as We increases, cj)+ decreases; here Pmif)) is the 
maximum value of P{4>). In particular, our DNS sug¬ 
gests that 1 — for small We. 

The modification in P{4>) can be understood qualita¬ 
tively by making the approximation that the effect of 
the fluid on the equation for (t> can be encapsulated into 
an eddy diffusivity Dg 
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The eddy-diffusivity- 
modified Cahn-Hilliard equation is dt4> = {De — D)V^(j) + 
DV^4>^ + MAW^cj), which gives the maximum and mini¬ 
mum values of 4> as (j)± = ^/ (D — Dp) /D. Furthermore, if 
we neglect the nonlinear term [l5l |20| , we find easily that 
the modified growth rate is Dk^[(l—De/D)—MAk^]\ i.e.. 


all wave numbers larger than kd = •\/(l — Dc/D)/{AM) 
are stable to perturbations. In particular, droplets with 
linear size <{2tt/ kd) decay in the presence of coupling 
with the velocity field; we expect, therefore, that, in the 
presence of fluid turbulence, the peak of P{4i) broadens 
and shifts as it does in our DNS. For a quantitative de¬ 
scription of this broadening and the shift of the peak, 
we must, of course, carry out a full DNS of the Cahn- 
Hilliard-Navier-Stokes equation as we have done here. 

We have investigated, so far, the effect of fluid turbu¬ 
lence on the phase-field (j) and its statistical properties 
such as those embodied in S{k) and P{(t>). We show next 


how the turbulence of the fluid is modified by (j), which 
is an active scalar insofar as it affects the velocity field. 
In the statistically steady state of our driven, dissipative 
system, the energy injection must be balanced by both 
viscous dissipation and dissipation that arises because of 
the interface, i.e., we must have Cinj = Ci/ + 

In Fig.|3|a), we show that decreases and increases 
as we increase We, while keeping Cinj constant, because 
Lc diminishes (Fig. [1]) and, therefore, the interfacial 
length and increase. This decrease of Lc is mirrored 
strikingly in plots of the fluid-kinetic-energy spectrum 
E{k) (Fig. mb)), which demonstrate that the inverse cas¬ 
cade of energy is effectively blocked at a wavenumber kc, 
which we determine below, from the energy flux, and 
which we find is ~ 2'k/Lc, where Lc follows from S{k) 
(see Fig. [2]). The value of kc increases with We; and the 
inverse cascade is completely blocked for the largest We 
we use, for which kc — kmj, the forcing scale. 

To provide clear evidence that the blocking of the 
energy flux is closely related to the arrest scale, we 
show in Fig. m c) plots of the energy flux Hsik) = 
T{k')dk' for different values of We. Here T{k) = 

Efc-i<fc'</c-Hi(«(-k,0-P(k) • (m X a;)(k,t))t is the en¬ 
ergy transfer and P(k) is the transverse projector with 
components Pij{k) = Sij — kikj/k"^. We define kc as 
the wave-number at which n£;(fc) becomes 4% of Cinj. 
We find that the wave-numer corresponding to the ar¬ 
rest scale 2tt/Lc (marked by vertical lines for each run) 
is comparable to kc. 


It has been suggested [4lj that coarsening arrest 
can be studied by using a model in which the field (j) is 
advected passively by the fluid velocity. Such a passive- 
advection model is clearly inadequate because it can¬ 
not lead to the phase-field-induced modifications in the 
statistical properties of the turbulent fluid (see Fig. El). 
However, for the sake of completeness, we now study the 
passive-advection case in which the coupling term (/A/g, 
is turned off in Eq. (ED. We then contrast the results for 
this case with the ones we have presented above. The 
parameters we use for the passive-advection DNS are 
N = 1024, A = = 0.0176; and we carry out runs 

for D = 5 - 10“^, 1 • 10“^, 5 • 10“^ and 5 • 10“^. The evolu¬ 
tion of the pseudo-grayscale plots of (f) with D, in the left 
panel of Fig. El is qualitatively similar to the evolution 
shown in Fig. [H There is also a qualitative similarity in 
the dependence on D of the scaled PDFs P{(j))/Pmic/); 
we can see this by comparing the passive-advection re¬ 
sult, shown in the middle panel of Fig. El for positive val¬ 
ues of in the vicinity of the peak, with its counterpart 
in Fig. El (c). However, there is a qualitative difference 
in the dependence of Lc on D: in the passive-advection 
case we find Lc ^ [Fig. El (inset)], which is in stark 

contrast to the essentially D-independent behavior of Lc 
shown in the inset of Fig. He). 


In conclusion, our extensive study of two-dimensional 






4 



FIG. 2. (Color online)(a) Log-log (base 10) plot of aLc versus e/cr'* showing data points {Lc from Eq. (|3]), with S{k) from our 
DNS) in red. The black line is the Hinze result Q for Lh', a fit to our data yields a constant of proportionality ~ 1.6 and an 
excellent approximation to the arrest scale Lc over several orders of magnitude on both vertical and horizontal axes; the plot 
of Lc versus D, in the inset, shows that, for fixed values of and cr (runs R1,R2 and R4), Lc is independent of D (within error 
bars), as is implied by the Hinze condition (see text), (b) Log-log (base 10) plots of the spectrum S{k), of the phase-field (j), 
versus fc; as VFe increases (i.e., cr decreases) the low-fc part of S{k) decreases and S{k) develops a broad and gentle maximum 
whose peak moves out to large values of k. (c) Plots versus (f), in the vicinity of the maximum at 0+, of the normalized PDFs 
P{cf))/Pm{4‘), where Pm{4>) is the maximum of P{4>); the peak position (j)+ —> 1 as We increases (see the inset which suggests 
that 1 — 0+ ~ (black line)). 


(2D) binary-fluid turbulence shows how the Cahn- 
Hilliard-Navier-Stokes coupling leads to an arrest of 
phase separation at a length scale Lc, which follows from 
S{k). We demonstrate that Lc ~ Lh, the Hinze scale 
that we find by balancing inertial and interfacial-tension 
forces, and that Lc is independent, within error bars, of 
the diffusivity D. We also elucidate how the coupling 
between the Cahn-Hilliard and Navier-Stokes equations 
modifies the properties of fluid turbulence in 2D. In par¬ 
ticular, we show that there is a blocking of the inverse 
energy cascade at a wavenumber kc, which we show is 
~ 27r/Lc. 


Earlier DNSs of turbulence-induced coarsening arrest 
in binary-fluid phase separation have concentrated on 
regimes in which there is a forward cascade of energy 
in 3D (see Ref. 23|) and a forward cascade of enstro- 
phy in 2D (see Ref. [22|). Although studies that use 
a passive-advection model for (j) obtain results that are 
qualitatively similar to those we obtain for S{k) and the 
spatiotemporal evolution of (j), they cannot capture the 
phase-field-induced modification of the statistical prop¬ 
erties of fluid turbulence and the correct dependence of 
Lc on D. We hnd our results to be in qualitative agree¬ 
ment with the earlier studies on advection of binary-fluid 
mixtures with synthetic chaotic flows ; of course, such 
studies cannot address the effect of the phase field on the 
turbulence in the binary fluid. 


Some groups have also studied the statistical properties 
of turbulent, symmetric, binary-fluid mixtures above the 
consolute point, where the two fluids mix even in absence 
of turbulence 4^ . In these studies, there is, of course, 

neither coarsening nor coarsening arrest. 

We hope our study will lead to new experimental stud¬ 
ies of turbulence in binary-fluid mixtures, especially in 
2D [lilil, to test the specific predictions we make for 


Lc and the blocking of the inverse cascade of energy. 

We thank S.S. Ray for discussions, the Department of 
Atomic Energy, the Department of Science and Technol¬ 
ogy, Council for Scientific and Industrial Research, and 
the University Grants Commision (India) for support. 
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line) with n_B(fc) gives fee, the wave-number at which the inverse energy cascade gets truncated; our estimate of arrest scale 
27r/Lc (vertical lines) is comparable to k^. 
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The parameters for our direct numerical simulations (DNSs) are given in Table HI 
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TABLE I: Parameters N, v, M, A, D, (fuiOj) for our DNS runs R1-R12. The forcing wave number is fixed at kinj = fiinj ~ 40, 
is the number of collocation points, v is the kinematic viscosity, M is the mobility, parameter, ^ sets the scale of the interface 
width, the surface tension cr = 2'\/2/3(A/^), (fujUi) is the enstrophy-injection rate, which is related to the energy-injection rate 
[cinj = ifuiO})/kf^j], E = O.SJDfel'W/sP is the fluid kinetic energy, is the fluid energy dissipation rate, 

€/_, = M is the energy dissipation rate because of the phase-field cf, is the coarsening-arrest length scale [Eq. (3)]. 

R1 — R2: Re = 124 and Sc = 1; R3 — RIO: Re = 124 and Sc = 0.1; Rll — R12: Re — 53 and Sc = 0.1. In some of our runs we 
also include a friction term —auj on the right-hand-side of Eq. (1); a = 0.001 for runs R4 — R8 and zero otherwise. The Cahn 
number Ch = ^/L = 






















